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Abstract. After a short review of electroweak scale baryogenesis, we consider the 
dynamics of chiral fermions coupled to a complex scalar field through the standard 
Yukawa interaction term at a strongly first order electroweak phase transition. By 
performing a systematic gradient expansion we can use this simple model to study 
electroweak scale baryogenesis. We show that the dominant sources for electroweak 
baryogenesis appear at linear order in the Planck constant K. We provide explicit 
expressions for the sources both in the flow term and in the collision term of the 
relevant kinetic Boltzmann equation. Finally, we indicate how the kinetic equation 
sources appear in the fluid transport equations used for baryogenesis calculations. 
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1 Introduction 

The necessary requirements on dynamical baryogenesis at an epoch of the early 
Universe are provided by the following Sakharov conditions: 

— baryon number (B) violation 

— charge (C) and charge-parity (CP) violation 

— departure from thermal and kinetic equilibrium 

The Sakharov conditions may be realised at the electroweak transition [1], pro- 
vided the transition is strongly first order. Namely, C and CP violation are re- 
alised in the standard model (SM) for example through the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix of quarks. B violation is mediated through the Adler- 
Bcll-Jackiw (ABJ) anomaly. At high temperatures the ABJ anomaly is mani- 
fest via the unsuppressed sphaleron transitions, and may be responsible for the 
observed baryon asymmetry today, which is usually expressed as the baryon-to- 
entropy ratio: 

^=3-7xlQ- u (1) 
s 

This is obtained both as a nucleosynthesis constraint and from recent cosmic 
microwave background observations. 

The standard model (SM) of elementary particles and interactions cannot 
alone be responsible for the observed matter-antimatter asymmetry (1), primar- 
ily because the LEP bound on the Higgs mass mg ^ 112 GeV is inconsistent 
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with the requirement that the transition be strongly first order. A strongly first 
order transition is namely required in order for the baryons produced in the 
symmetric phase not be washed-out by the sphaleron transitions in the Higgs 
('broken') phase. And this is so provided the transition is strong enough. This 
is usually expressed as the requirement A<p J; 1.1T on the jump in the Higgs 
expectation value <f> of the phase transition [2]. 

Supersymmetric extensions of the Standard Model on the other hand may 
result in a strongly first order transition. For example, in the Minimal super- 
symmetric standard model the sphaleron bound can be satisfied provided the 
stop and the lightest Higgs particles are not too heavy, ~ 1 20 GeV and 
m H <, 120 GeV [3]. 

An efficient mechanism for baryon production at the electroweak phase tran- 
sition is the charge transport mechanism [4], which works as follows. At a first 
order transition, when the Universe supercools, the bubbles of the Higgs phase 
nucleate and grow. In presence of a CP-violating condensate at the bubble in- 
terface, as a consequence of collisions of chiral fermions with scalar particles in 
presence of a scalar field condensate, CP- violating currents are created and trans- 
ported into the symmetric phase, where they bias baryon number production. 
The baryons thus produced are transported back into the Higgs phase where 
they are frozen-in. The main unsolved problem of electroweak baryogenesis is 
systematic computation of the relevant CP-violating currents generated at the 
bubble interface. Here we shall reformulate this problem in terms of calculating 
CP-violating sources in the kinetic Boltzmann equations for fermions. 

tree =—^-t — + —> — ;— + . : ■ + .. 
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x CP violating condensate 

Fig. 1. The tree level interactions of fermions with the scalar field condensate which, 
when expanded in gradients, lead to the CP- violating semiclassical force in the kinetic 
Boltzmann equation. 

The techniques we report here are relevant for calculation of sources in the 
limit of thick phase boundaries and a weak coupling to the Higgs condensate. 
In this case one can show that, to linear order in the Planck constant h~, the 
quasiparticle picture for fermions survives [5, 6]. In presence of a CP-violating 
condensate there are two types of sources: the semiclassical force in the flow term 
of the kinetic Boltzmann equation, and the collisional sources. The semiclassi- 
cal force was originally introduced for baryogenesis in two-Higgs doublet models 
in [7], and subsequently adapted to the chargino baryogenesis in the Minimal 
Supersymmetric Standard Model (MSSM) in [8]. The semiclassical force cor- 
responds to tree level interactions with the condensate shown in figure 1 and 
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it is universal in that its form is independent on interactions. The collisional 
sources on the other hand arise when fermions in the loop diagrams interact 
with scalar condensates. In figure 2 we show typical CP-violating one-loop con- 
tributions to the collisional source. This source arises from one-loop diagrams in 
which fermions interact with a CP-violating scalar condensate. When viewed in 
the kinetic Boltzmann equation, these processes correspond to tree-level inter- 
actions in which fermions absorb or emit scalar particles, whilst interacting in 
a CP-violating manner with the scalar condensate. The precise form of the col- 
lisional source depends on the form of the interaction. In the following sections 
we discuss how one can study the CP-violating collisional sources induced by a 
typical Yukawa interaction term. 
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Fig. 2. The one- loop fermion-scalar diagrams, where interactions of fermions with the 
scalar field condensate are explicitly shown. The condensate interactions, when ex- 
panded in gradients, result in a CP- violating collisional source in the kinetic Boltzmann 
equation. Upon projecting (or 'cutting') the loop propagators on-shell one obtains the 
CP-violating scalar particle absorption and emission processes. 



2 Kinetic equations 

Here we work in the simple model of chiral fermions coupled to a complex scalar 
field via the Yukawa interaction with the Lagrangian of the form [5, 6] 

C = i^fiip - Tp L mip R - ip R m*Tp L + £ Yu , (2) 

where £v u denotes the Yukawa interaction term 

£y u = -vH^LipR - y<t>*i>Ri>L, (3) 

and m is a complex, spatially varying mass term 

m(it) = y ,( Po = n^niu) + imi(u) — \m(u)\e lS ^ u \ (4) 

Such a mass term arises naturally from an interaction with a scalar field con- 
densate <Po = (@(u)). This situation is realised for example by the Higgs field 
condensate of a first order electroweak phase transition in supersymmetric mod- 
els. When <t> m (3) is the Higgs field the coupling constants y and y' coincide. 
Our considerations are however not limited to this case. 
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The dynamics of quantum fields can be studied by considering the equations 
of motion arising from the two-particle irreducible (2PI) effective action [9] in 
the Schwinger-Keldysh closed-time-path formalism [10, 11]. This formalism is 
suitable for studying the dynamics of the non-equilibrium fermionic and bosonic 
two-point functions 

iS a0 (u,v) = {Q\T c [Mu)M v )\ W) ( 5 ) 
iA(u,v) = (f2\T c [<j>(u) ( f>Hv)]\{2), (6) 

where \Q) is the physical state, and the time ordering Tq is along the Schwinger 
contour shown in figure 3. For our purposes it suffices to consider the limit when 
to — > — oo. The complex path time ordering can be conveniently represented in 
the Keldysh component formalism. For example, for noncquilibrium dynamics 
of quantum fields the following Wightman propagator is relevant 

iS < (u, v) = ~{Q\ip(v)ip(u)\n). (7) 



% C+ 

Fig. 3. The complex time contour for the Schwinger-Keldysh nonequilibrium formal- 
ism. 



For thick walls, that is for the plasma excitations whose de Broglie wavelength 
£dB is small in comparison to the phase interface thickness L w , it is suitable to 
work in the Wigner representation for the propagators, which corresponds to the 
Fourier transform with respect to the relative coordinate r = u — v, and expand 
in the gradients of average coordinate x = (u + v)/2. This then represents an 
expansion in powers of £d&/L w . When written in this Wigner representation, the 
kinetic equations for fermions become [12] 

VS< = (i;9 + #- (mP R -m*P L )e-^ d d ^S < = (8) 

where for simplicity we neglected the contributions from self-energy corrections 
to the mass and the collisional broadening term [12]. When the collision term 
C$ is approximated at one loop, equation (8) corresponds to the noncquilibrium 
fermionic Schwinger-Dyson equation shown in figure 4. Since the scalar equation 
(also shown in figure 4) does not yield CP-violating sources at first order in 
gradients [5, 14], we shall not discuss it here. 

As the bubbles grow large, they tend to become more and more planar. 
Hence, it suffices to consider the limit of a planar phase interface, in which the 
mass condensate in the wall frame becomes a function of one coordinate only, 
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Fig. 4. The one-loop Schwinger-Dyson equations for the out-of-equilibrium fermionic 
(S) and scalar (A) propagators. When projected on-shell and expanded in gradients, 
these equations reduce to the kinetic Boltzmann equations. 

m = m{z). Further, we keep only the terms that contribute at order h to Eq. (8), 
which implies that we need to keep second order gradients of the mass term 

me -i a- d k = m + Im'dk, - \m"dl + o(fl£), (9) 
2 8 

where m — m(z), m! = d z m and m" = d z m. On the other hand, in the collision 
term Q, we need to consider terms only up to linear order in derivatives 

C^p = C,po + C^,i + .. 

C0i = - dj»d?))(E>S< E<S>), (10) 

where S < and A7 > represent the fermionic self-energies, and the derivatives 
@k (^ 2 \ djp) act on the first (second) factor in the parentheses. 

An important observation is that, when G = Gik^, t — xn • ky , z), the spin in 
the z-direction (of the interface motion) 

S z = L- 1 (A)S Z L(A) = 7|| (§ t - i( V || x a),) (11) 

is conserved 

[X>,S*]S< =0, (12) 

where £> is the differential operator in Eq. (8), a — 7°7, S z — 7°7 3 7 5 , and 
7|| = 1/(1 — Vn) 1 / 2 . This then implies that, without a loss of generality, the 
fermionic Wigner function can be written in the following block-diagonal form 

s=± 

S< =L(A)- 1 S < L(A) 
~i7°S< = ±{l + sa 3 )®p a g s a , (13) 



6 Tomislav Prokopec 



where <r 3 and p 1 (i — 1, 2, 3) are the Pauli matrices and p° = 1 is the 2x2 unity 
matrix and L(A) is the following Lorentz boost operator 

L(A) = fc0 + ^ Vk ", (14) 
^/2fc (fco + fco) 

with fc = sign(/c )(fco — l* 2 ) 1 / 2 - The boost A corresponds to a Lorentz transfor- 
mation that transforms away kii . 

With the decomposition (13) the trace of the antihcrmitean part of Eq. (8) 
can be written as the following algebraic constraint equation [6] 




where gQ — 7||ffo < denotes the particle density on phase space {k^^Xv}. Equa- 
tion (15) has a spectral solution 

•9oo = ^Z^~ n s5{k a tw,±), (16) 

_|_ -^s± 

where ui s ± denotes the dispersion relation 



w s ±=w ts- — — , uj = Jk 2 + H 2 , wo = Jul- k, 2 (17) 

Zujqujq v 11 

andZ s ± = l^s\m\ 2 9' /2u)q. The delta functions in (16) project n s (fc M , i — X|| -k|| , z) 
on-shell, thus yielding the distribution functions f s+ and / s _ for particles and 
antiparticles with spin s, respectively, defined by 

/ s+ =n s (w s+ ,fc z ,i-X|| -k||,z) 

/ s _ = 1 - n s (-u) s -, -k z ,t + X|| • k||,z). (18) 

This on-shcll projection proves the implicit assumption underlying the semi- 
classical WKB-mcthods, that the plasma can be described as a collection of 
single-particle excitations with a nontrivial space-dependent dispersion relation. 
In fact, the decomposition (13), Eq. (15) and the subsequent discussion imply 
that the physical states that correspond to the quasiparticle plasma excitations 
are the eigenstates of the spin operator (11). 

Taking the trace of the Hermitean part of Eq. (8), integrating over the pos- 
itive and negative frequencies and taking account of (16) and (18), one obtains 
the following on-shell kinetic equations [6] 

dtfs±+v\\ ■ V||/ s ± +v s ±d z f s ± +F s ±d k J s ± =C i , s ±[f s ±], (19) 

where f s ± = f s ±{k ll ,z, t — V| • X||), C s ±[f s ±] is the collision term obtained 
by integrating (10) over the positive and negative frequencies, respectively, the 
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quasiparticlc group velocity v s ± = k z /ui s ± is expressed in terms of the kinetic 
momentum k z and the quasiparticle energy w s ± (17), and the semiclassical force 

F^-tf ± fW. (20) 

ZOJ s ± ZLOqLOq 

In the stationary limit in the wall frame the distribution function simplifies to 
f s ± — f s ±{kfj,, z, ). When compared with the 1+1 dimensional case studied in [5] 
the sole, but significant, difference in the force (20) is that the CP-violating 9'- 
term is enhanced by the boost-factor 7|| = wq/u)q, usq = {oj^ — k^ 2 ) 1 / 2 , which, 
when integrated over the momenta, leads to an enhancement by about a factor 
two in the CP-violating source from the semiclassical force. 



3 Sources for baryogenesis in the fluid equations 

Fluid transport equations are usually obtained by taking first two moments 
of the Boltzmann transport equation (19). That is, integrating Eq. (19) over 
the spatial momenta results in the continuity equation for the vector current, 
while multiplying by the velocity and integrating over the momenta yields the 
Eulcr equation. The physical content of these equations can be summarized as 
the particle number and fluid momentum density conservation laws for fluids, 
respectively. This procedure is necessarily approximate simply because the fluid 
equations describe only very roughly the rich momentum dependence described 
by the distribution functions of the Boltzmann equation (20). The fluid equations 
can be easily reduced to the diffusion equation which has so far being used almost 
exclusively for electroweak baryogenesis calculations at a first order electroweak 
phase transition. A useful intermediate step in derivation of the fluid equations 
is rewriting Eq. (19) for the CP-violating departure from equilibrium 5f si = 
Sf si+ - 5f si - as follows 



d t + —d z 



dk z Sf si + v w 5F si (d u f u ) Uoi 



+ v w F iSw s 



Fqi — - 
Su) si = s 



TO,; 



■"ipsit 



^uj fid 

where i is the species (flavour) index, f u = l/(e^ w + 1), and 

12' 



$Fsi — Fsi-\- -^si— 5 



■"4>si 



ipsi- 



(21) 



(22) 
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When integrating (21) over the momenta, the flow term yields two sources in 
the continuity equation for the vector current. The former comes from the CP- 
violating spin dependent semiclassical force, and has the form 



-sv u 



d 3 k 
(2tt) 3 

(\mi\XY 



4tt 2 



3-F s i {duj fuJ^LU — LJQi 

\J a \p^i ) 



(23) 



with Xi — \rm\/T, while the latter comes from the CP-violating shift in the 
quasiparticle energy, and can be written as 



— 



d 3 k 

W) 



:Foi8ci! s i 



(|TOi| 2 6>-) , 

' vfr- \mi\ Jb(xi). 



(24) 



The total source is simply the sum of the two, S S i = + . To get a more 




1 2 3 4 5 6 



Xj 

Fig. 5. The flow term sources (23)-(24) characterised by the integrals xlJ a (xi) {red 
solid) and x\jh(xi) {green dashed) as a function of the rescaled mass Xi = \m,i\/T. The 
sum of the two sources {thin blue line) is also shown. 



quantitative understanding of these sources, in figure 5 we plot the integrals J a 
and Jb in equations (23) and (24). A closer inspection of the sources <S^ and 
S h si indicates that the total source S S i can be also rewritten as the sum of two 
sources: the source oc Im^ 2 '©-, characterized by x\j a + x 3 Jb, and the source 
oc |mi| 2 (9", characterized by xfj a . We note that in the spin state quasiparticle 
basis the flow term sources appear in the continuity equation for the vector 
current, while in the helicity basis, which is usually used in literature [7, 8], the 
flow term sources appear in the Euler equation. In figure 6 we show recent results 
of baryogenesis calculations of Ref. [8] based on the CP-violating contribution 
to the semiclassical force in the chargino sector of the Minimal Supersymmetric 
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Standard Model (MSSM) . This calculation is based on the quasiparticle picture 
based on helicity states. The analysis suggests that one can dynamically obtain 
baryon production marginally consistent with the observed value (1), provided 
m-2 ~ (X ~ 150 GeV and v w ~ 0.03c. 
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Fig. 6. The semiclassical force baryogenesis mediated by charginos of the MSSM cal- 
culated in the helicity basis. The figures show contours for the baryon-to-entropy ratio 
in the units of 10 -11 for two wall velocities v w = 0.01 and v w = 0.03 as a function of 
the soft susy breaking parameters fi and mi. A maximal CP violation in the chargino 
sector is assumed. The shaded (yellow) regions are ruled out by the LEP measurements. 
The observed baryon asymmetry is in these units 3 — 7. (The figure is taken from the 
latter reference in [8]. ) 



We now turn to discussion of the collision term sources in Eqs. (19) and (21). 
We assume that the self-energies S >,< are approximated by the one-loop expres- 
sions (cf. figure 4) 

E<>{k, x) = iy 2 J d4 ^" [{2*f8{k — k' + k")P L S<>>{k', x)P R A>-<{k",x) 
+ (2ir) 4 5(k - k' - k")P R S < > > {k',x)P L A < ' > (k",x)] , (25) 

where Z\ < and Z\ > denote the bosonic Wigner functions. This expression con- 
tains both the CP-violating sources and relaxation towards equilibrium. The 
CP-violating sources can be evaluated by approximating the Wigner functions 
S > ' < and A >,K by the equilibrium expressions accurate to first order in deriva- 
tives. The result of the investigation is as follows. There is no source contributing 
to the continuity equation, while the source arising in the Euler equation is of 
the form [12] 

f d 4 k k z 2 S H 261 ' 
2 J ± j2^- C ^ = VwV 32W^ (|m| ' m ^' (26) 

where the function Tf(\m\, m^) is plotted in figure 7. It is encouraging that the 
source vanishes for small values of the mass parameters, which suggests that the 
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expansion in gradients we used here may yield the dominant sources. Note that 
the source is nonvanishing only in the kinematically allowed region, > 2\m\. 
When the masses are large, \m\,m < f > 3> T, the source is as expected Boltzmann 
suppressed. It would be of interest to make a detailed comparison between the 
sources in the flow term and those in the collision term. This is a subject of an 
upcoming publication. 
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Fig. 7. The collisional source contributing to the fermionic kinetic equation at one loop 
for the mass ratios m^/lml = 2.1,2.5,3,4,10 and 20, respectively. The source peaks 
for \m\ « 0.7T and w 4|m|. 
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